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I. Summary of Research Under NASA Grant NGR14-005-009 
1*1. Brief Review 

The research on the NASA Grant grew out of a rocket experiment on 
the cross-modulation between two RF signals in the ionosphere, the 
so-called "Luxemburg effect." In the original experiment over ten years 
ago, a pulsed signal near the gyroresonance frequency was produced to 
heat the electrons in the E layer, and a sensing signal was sent through 
the heated region for detection of the heating effect. The local dis- 
turbance of the ionosphere was observed. However, it was found that an 
accurate quantitative interpretation of the experimental results required 
a thorough understanding of the behavior of antennas in anisotropic media. 
The problem of the antenna in magnetoplasma has been a subject of great 
interest in the past decade because of its intimate relevance to space 
communication. During the past few years our laboratory has contributed 
significantly in this area as shown in the list of publications in 1.2. 
Some of the highlights of this research is now reviewed. 

(1) A theory for a short dipole antenna in a magnetoplasma was 

4 * 

advanced and a formula for the antenna impedance was 

derived using the quasi-static approximation. This formula 

has been widely used and experimentally tested with some 

4 36 

success both in the ionosphere and in the laboratory. * 

(2) The computation of Green f s function in the magnetoplasma 
has been a problem because of its complexity. A significant 
simplification was achieved with a new expression which 
separates the singular algebraic terms from some proper 
integrals. 1 More recently this expression has been 
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rearranged so that it can be handled more easily by some 


» 


numerical means (see (6) below) . 

(3) For arbitrary electric current sources, the far field can be 

expressed in terms of the geometrical properties of the 

medium dispersion surface."* 5 ^ These results were later 

15 c 

extended to compressible magnetoplasma and to sources or 
magnetic currents, mechanical forces, and actual injection 
of charged particles. The possible significance of these 
sources in the case of compressible plasma was demonstrated, 
using a variational method. 

(4) Rigorous and numerical solutions to the current distribution 

on an infinitely long cylindrical antenna parallel to the 

23 

magnetizing field were obtained for various cases. These 


17,25 


solutions differ from others in that the boundary condition 
on the conducting cylinder was satisfied. 


(5) Dipole impedance in a plasma-filled circular guide has been 
investigated both numerically and experimentally with good 
agreement. 24 It was shown that the interaction between the 
dipole and the plasma could lead to a useful diagnostic 
technique for the plasma. 

(6) In the past the investigations were made largely on the 
assumptions that the medium was idealized; the sheath was 
neglected or represented by a dielectric layer, and the current 
distribution was known (except (4) where the antenna was 
assumed to be infinitely long) . Xn the past year we have 
addressed ourselves to the problem of current distribution. 
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The method relies heavily on the computer. It would have 
been unthinkable a few years ago, but with the present 
availability of larger and faster computers we have found 
this approach feasible with a reasonable computing effort. 

An integral equation was solved numerically for the 
antenna current distribution in some cases. Because of the 
extreme complexity of the Green’s function, the numerical 
procedure is very involved. The program has been successfully 
tested for simple known cases. 

(7) In the past year, an experimental program was also initiated 
for the study of the nonlinear effect of the plasma on the 
antenna. It was found that even for input power on the order 
of milliwatts there was significant nonlinear effect on the 
antenna impedance. This should be kept in mind in evaluating 
the results of many existing experimental works. The nonlinear 
effect of the plasma depends on its relaxation time in comparison 
with the periods corresponding to the signal and modulation 
frequencies. At microwave frequencies, it is an extremely 
complex problem, and rigorously speaking the antenna impedance 
as well as the so-called ei-characteris tics in such a system 
become rather meaningless. Therefore a more complete description 
is needed. Since at microwave frequencies it is usually the 
measurable reflection coefficient that is of interest, its 
dependence on the level of incident wave in the case of a 
simple nonlinear plasma has been evaluated. 

More recently other nonlinear effects, such as harmonic 
generation and cross— modulation between two or more signals, 
have been demonstrated in the laboratory plasma. (See Part II 


of the report) . 
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I. 2. List of Publications and Scientific Reports under NASA Grant 
NSG 395 and NGR 14-005-009 


1 . 


2 . 


3 . 


4 . 


5 . 


6 . 


7 . 



9 . 


10 . 


11 . 


12 . 


13. 


R. Mittra and G. A. Deschamps, "Field Solution for a 
Anisotropic Medium, " Proceedings of the Internatioiial 
Electromagnetic Theory, Pergamon Press, pp. 495-512, 


Dipole in an 
Symposium on 


1963. 


G. A. Deschamps. "Dispersion Surfaces and Characteristic Plane Waves 
in a Gyrotropic Medium," Proceedings of the Application Forum on 
Antenna Research, ed. P. E. Mayes, University of Illinois, pp. 290-367, 
January 1964. 

R. Mittra, "Antennas in Anisotropic Media," Applications Forum on 
Antenna Research, ed. P. E. Mayes, pp. 424-452, January 1964. 

K. G. Balmain, "The Impedance of a Short Dipole Antenna in a Magneto- 
plasma," IEEE Transactions, Vol. AP-12, pp. 605-617, September 1964. 


G. A. Deschamps and 0. B. Kesler, "Radiation Field of an Arbitrary 
Antenna in a Magnet op 1 as ma," IEEE Transactions, Vol. AP-12, pp. 7S3- 
785, November 1964. 


G. A. Deschamps, "Angular Dependence of the Refractive Index in the 
Ionosphere," Radio Science, 69D, pp. 395-400, March 1965. 

K. G. Balmain, "Impedance of a Short Dipole in a Compress iblc Plasma," 
Radio Science, 69D, No. 4, pp. 559-566, April 1965. 

R. Mittra and'G. L. Duff, "A Systematic Study of the Radiation Patterns 
of a Dipole in a Magnetopiasma Based on a Classification of the 
Associated Dispersion Surfaces," Radio Science, 69D, No. 5, pp. 631-69a, 
May 1965. 

C. P. Bates and R. Mittra, "Plane-Wave Spectrum Representation of a 
Field in a Magnetoionic Medium Due to a Plane Distribution of Source 
Currents," Electronics Letters', Vol. 1, No. 3, p. 74, May 1965. 


S. W. Lee, C. Liang, and Y. T. Lo , "Inconsistency of Boundary Conditions 
of Plasma Models in a Bounded Region," Electronics Letters of the IEEE, 
Vol. 1, No. 5, p. 128, July 1965. 


S. W. Lee and Y. T. Lo, "Radiation Resistance of an Elliptical Loop 
Antenna with Constant Current in Compressible Plasma," Electronics 
Letters of the IEEE, Vol. i. No. 5, pp. 132-133, July 1965. 


G. L. Duff and E. Mittra, "Input Impedance of a Small Loop of Uniform 
Electric Current in an Anisotropic Cold Plasma," Electronics Letters, 

Vol. 1, No. 5, pp. 117-118, July 1965. 

S. W. Lee, C. Liang and Y. T. Lo, "Further Remarks on Boundary Conditions 
of Plasma Models in a Bounded Region," Electronics Letters of the IEEE, 
Vol. 1, No. 7, September 1965. 
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14. S.1V. Lee, Y.T. Lo and R. Mittra, "Finite and Infinite H-plane Di- 
furcation Waveguide with Anisotropic Plasma Medium," Proceedings of 

the International Symposia of Antenna and Propagation , Washington/ D. C. 
August 196S; also Canadian Journal of Physics , Vol. 43, pp. 2123-2135, 
December 1965. 

15. G. A. Desehamps, "Rayonnement d’une Antenne dans un milieu anisostrope," 
L’Onde Electrique 45 , pp. 1379-13S5, December 1965. 

16. S. W. Lee and Y. T . Lo, "Radiation in a Moving Anisotropic Medium," 
Radio Science Vol. 1, (New Series), No. 3, pp. 313-324, Marcii 1966. 

17. J. W. Carlin and R. Mittra, "Terminal Admittance of a Thin Biconical 
Antenna in an Isotropic Compressible Plasma," Electronics Letters , 

Vol. 2, No. 6, pp. 199-201, June 1966. 

18.. S. W. Lee, C. Liang and Y. T. Lo, "Wave Propagation in Plasma with Very 
Strong Magnetic Field," Radio Science , Vol. 1 (New Series), No. 7, 
pp. 815-824, July 1966. 

19. C. P. Bates and R. Mittra, "Principles of Scaling in a Uniaxial Medium," 
II Nuovo CimentOj Vol. 43, No. 10, pp. 101-118, 1966. 

20. G. A. Desehamps, "Impedance of Antennas in Ionized Media," Electro- 
magnetics and Antennas > ed. P. E. Mayes, University of Illinois, 
pp. 498-512, January 1967. 

21. S. W. Lee, "Propagation in Ionosphere-earth Waveguide," Appl. Science 
' Research , Vol. 17, No. 2, pp. 115-132, 1967. 

22. S. V/. Lee and Y. T. Lo, "Reflection and Transmission of Electromagnetic 
Waves by a Moving Anisotropic Medium," Journal of Applied Physics , 

Vol. 38, No. 2, pp. 870-875, February 1967. 

23. S. W. Lee and Y. T. Lo, "Current Distribution and Input Admittance of 
a Cylindrical Antenna in Anisotropic Plasma," presented at Spring URSI 
Meeting, Washington, D. C., April 1966; also, IEEE Transactions > Vol. 
AP-15, No. 2, pp. 244-252, March 1967. 

24. S. W. Lee and Y. T. Lo, "On the Coupling of Modal Waves in a Plasma- 
filled Parallel-plate Waveguide," Radio Science j Vol. 2, No. 4, 

pp. 401-406, April 1967. 

25. J. W. Carlin and R. Mittra, "Acoustic Waves and Their Effects on 
Antenna Impedance," IEEE Inti . Antennas and Propagation Symposium] 

Digest , pp. 102-105, December 1966, also, Canadian Journal of Physics > 
Vol. 45, pp. 1251-1269, 1967. — ; 

26. S. W. Lee and R. Mittra, "Group Velocity of a Uniaxial Plasma," IEEE 
Transactions^ Vol. AP-15, No. 4, pp. 588-589, July 1967. 

27. G. A. Desehamps and 0. B. Kcsler, "Radiation of an Antenna in a Com- 
pressible Magnetoplasma," Radio Science^ August 1967. 
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28. R. Mittra and S. IV. Lee, "Mode Matching Method for Anisotropic Guidos,’- 1 
Radio Science 3 Vol . 2, No. 8, pp. 937-942, 1967. 

29. S. W. Lee and R. Mittra, "Transient Radiation of an Electric Dipole 
in a Uniaxially Anisotropic Plasma," Radio Science } Vol. 2, No. 3, 
pp. 813-820, 1967. 

30. R. Mittra and S. IV. Lee, "Discontinuity in an Anisotropic Waveguide," 
Journal of Applied Physics * Vol. 38, pp. 178-186, 1967. 

31. S. IV. Lee a*nd R. Mittra, "Edge Coiidition and 'Intrinsic Loss’ in 
Uniaxial Plasma," Canadian Journal of Physics , Vol. 46, pp. 111-120, 
1968. 

32. J. Carlin and R. Mittra, "Effects of Induced Acoustic Sources on the 
Impedance of a Cylindrical Dipole in a Warm Plasma," Radio Science , 

Vol. 2, No. 11, pp. 1327-1538, 1967. 

33. J. Carlin, "Variational Formula for Antenna Impedance in a Warm Plasma," 
Electronics Letters , Vol. 3, No. 7, 1967. 

34. D. K. Waineo and R. Mittra, "Short. Dipole in a Uniaxial Medium," Applied 
Science Research , Vol. 18, 1968. 

35. I. Akkaya, "A Comparative Study of Several Models of the Magnetoplasma," 
Radio Science > Voi. 4, No. 1, pp. 83-90, January 1969. 

36. D. Snyder and R. Mittra, "Experimental Study of the Impedance of a 
Short Dipole in a Plasma for Parallel and Perpendicular Orientation 
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Science , Vol. 4, No. 2, pp. 179-189, February 1969. 
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39. G. L. Duff and R. Mittra, "Loop Impedance in a Magnetoplasma: Theory 

and Experiment," Radio Science 3 Vol. 5, pp. 81-94, January 1970. 

^0* T. I toh and R. Mittra, "A New Method of Solution for Radiation from 
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II. 1. Summary 

The research during the last period consisted of computations 
toward the solution of the problem of the current distribution on a 
cylindrical antenna in a magneto plasma. The case of an antenna parallel 

I 

to the applied magnetic field was treated (reference 63 of Part I) . A 
systematic method of asymptotic expansion was found which simplifies 
the solution in the general case by giving the field of a dipole even 
at relatively short range (II. 2). Some useful properties of the dis- 
persion surfaces in a lossy medium have also been found (II. 2). The 
laboratory experiment was directed toward evaluating nonlinear effects, 
such as those due to power level, bias voltage and electron heating (II. 3). 
The problem of reflection and transmission of waves in an electron heated 
plasma was treated theoretically (II. 4). [This has been accepted for 

publication]. The profile inversion problem has been pursued. Some 
results are very encouraging. The general question of stability of the 
solution remains unsolved. 
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II. 2 . Cylindrical Dipole in a Magnetoplasma (Numerical Analysis) 

1. Introduction 

The input impedance to a dipole in a magnetoplasma has been computed 
using the method of moments. This paper discusses asymptotic approxi- 

mations to integrals representation of fields. By applying these approxi- 
mations the computation of fields from a dipole source in a plasma can be 
simplified and still give good results at "reasonable" distance from the 
source. In this manner, the asymptotic solution can b.e. applied to the 
problem of obtaining the input impedance of a dipole in a magnetoplasma. 
Such an application could increase the accuracy of the computation. 

This, paper discusses methods for obtaining the first few terras of 
an asymptotic expansion. Both one and two dimensional integrals are 
considered of the type 

/ (weighting function) (plane wave) d (parameter) (1) 

c 

The field solution in the free space case is examined as a preliminary 
step to the eventual application to a plasma. Also some geometrical 
considerations about the characteristic fields in an anisotropic region 
are presented. 

2. Saddle Point Integration 

r 2 ] 

The saddle-point method of integration, also known as the method 
of steepest descent, is an asymptotic integration technique which may have 
some usefulness in integrating the equation 

1 - f f (x) e 1V( ^ x) dx (2) 
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where v is a large number and f, (j> are functions of x that vary slowly 
with x ■— x being a point in a space of arbitrary dimension. The particular 
application of ( 2 ) used here is the consideration of the field described 
by a spectrum of plane waves. Then the variable x is the vector k, the 
phase V$(x) becomes k * r (where r is the observation point) , and F(x) 
becomes a weighing function F(k). The integral considered thereby becomes 

/ F(k)e lk r dk (3) 

c 

The path of integration c is shifted to a steepest descent path that is 

one along which Im(k • r) varies at a maximum rate of change. This can 

be done without changing the value of the integral only if there are no 

singularities between the true path. Otherwise the effect of the poles 

or branch points on the new integration path must be considered. 

r 3 1 

It is easy to show 1 J that forcing Im(k • r) to vary at a maximum 
rate. of change implies that Re(k * r) is a constant. 

A complex analytic function cannot have a maximum or minimum } so 
the stationary phase point is actually a saddle point. There is both 
a steepest descent path, along which Re(k • r) has its maximum at the 
saddle point, and also a steepest ascent path along which Re (k • r) 
has its minimum at the saddle point. 

If one can change the integration path of (3) to the steepest descent 
path and if F(k) varies "slowly” one now has a means for approximating 
the solution to (3). Specifically, if one integrates along the steepest 
descent path then the main contribution to the integration comes from the 
region of the saddle point. In this manner expanding the integral in 
terms of a Taylor’s series about the saddle point for "well-behaved" 
functions should produce reasonable approximations to the solution of (3) . 


r~ 



13 


At a saddle point d(kr) = 0 or for a given r, dk • r = 0. This 
means that r is perpendicular to the dispersion surface. At point k, 
in trying to find an approximation to the integral solutions to a field 
at an observation point r one must find the saddle points. This intro- 
duces a geometric problem necessary for the solution of the integral (3) 
of finding all those k's for which dk • r = 0. This problem will be 
discussed later. 

Consider now a one-dimensional problem with an integral of the form 


I - / F( Z ) e " af(z> d Z 


'sdp 


w 


Here is the steepest descent path. It is assumed that F(z) varies 

"slowly" and a is "reasonably" large. Also one assumes that F(z) is 
well-behaved at the saddle-point, having no singularities there. The 
approximation to (4) about one saddle point is now examined. Two methods 


are mentioned here. 

2 

In approximation (4) define s - f - f(z^) where z^ is the saddle 
point. Making a change of variables from z to s results in 


-a f(z ) 
1 - e o 


$(s) 


2F(z) 

f'W 


/ $(s) 

C sdp 



(5) 


To evaluate (5), $(s) is expanded about the point s = 0 (or z = Z Q )* 
th 

Let the n derivative of $(s) and f(s) evaluated at the saddle point 
be denoted by and f respectively. The Taylor’s expansion for <{)(s) 
becomes 


2 ‘ 

<f>(s) = 4 > q + + ^ 2 * ^2 + 


( 6 ) 
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and <j> = 0 at the saddle point. In order to have the Im (-as ) a constant, 

2 

s must be real. So, appropriate approximations to the steepest descent 
path valid about the saddle point for <f>(s) varying "slowly” are limits 
of integration - 00 to +°°. Putting the expansion (6) into (5) yields an 
approximate solution to (5) of 


I = / I e af o [<J, o + ^ <f> 2 + * 4 ...] 

8a 


(7) 


This computation has been performed to the first two significant terms 


using powers of i and F n with the following results 


I = e 


-af ,2 tt x ‘ L/iL ' / ^ 2 "h 5 "3 N , 1 *3 „ ^ 

0 fe 0 {F ‘ 37 ( f7 ' 3 3 > + 4 (F 1 £7 ' V } 


1/2 


f 

2 f 2 


. -af , ,2 , / f 4 x 2 35 

+ e o (- ^af ^ 4f 3 48 ( f^ 1536 


2 

f 3 f 5 7 f 3 f 4 7 f 3 £ 4 35 

2 192 2 192 3 256 

t 2 1 2 t 2 


+ F rfi 1 .. V4_35 + ' f f3,3_35 

1 l f_ 32 ,2 192 v f ' 192 J 

2 f „ 2 


F f 

+ 4t 4 5 


385 


f 3 4 

(■A) 


2 l f 2 32 (64) (72) v f 2 


f 2 

(- 1 ) 35 ] 

f 2 96 


- F 3 Is + F 4 32 3 } 


( 8 ) 


An equivalent expansion for the one-dimensional problem has been 
found by a slightly different method which is more readily extended to 
the nth dimensional case. In equation (4), F(z) and f(z) are expanded 
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in a Taylor’s series about a saddle point z. Let £ - z - z q . Then 


I = / (F q + v ± C + -f C 2 + ...) 
C sdp 


(9) 


exp 


2 J2 


(-a{f + f, ? + -4 C + ...» d? 

o i z 


The exponential term in (9) is now expanded. 

f 

exp (-a f(z)) = exp(-af) exp(- Z?) 

^3 3 ^4 4 ^5 ^5 

{1 - a(— C + -ry C + r; • * •) 


5! 


2 f ~ « f/ . fr c 2 

+ il_ (--! +■ — . , , ) 

^3! ^ 4! ^ 5! ^ J 


2! v 3 ! ■" ' 4! 

3 3 

ff (....) ...)} 


( 10 ) 


substituting (10) into (9) reduces the integral equation to an infinite 
sum of terms of the form 


-af , c r n e~ ^ ^ 

e o / c n ? e ~ 

C sdp 


d e 


(ii) 


where C is a constant. The limits of integration are taken from -«> to 00 
n 

as an approximation around the saddle point. 

With this simplification, it is known that 

r r[(n + l)/2] . 


» 2 2 
/ C" e"P « d C = 


(n + l)/2 


n even 


n odd 


( 12 ) 


1 1 ^ 2 1 

When a sum of these factors is examined in terms of — and — , a 

a a 


series equivalent to (5) has been obtained. While this method is more 
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cumbersome because of the number of terms that must be combined, it is 
a straight forward technique that may easily be applied to the two- 
dimensional problem. 


3. Two-Dimensional Problem 


In a two-dimensional problem, one considers integrals of the type 


I 2 = ff F(u,v) e " pf(u ' v) du dv (13) 

C sdp 

Again it is assumed that F(x,y) varies "slowly" and is well-behaved at 
the saddle point. 

To simplify the manipulation of (13) a convenient notation is first 

introduced. Let -J be an ordered pair of indices so that J “ ^1^2' 

2 . 

Also, let y be a vector in C so that y^ - u and y^ = V. The following 

nbtation will now be used. 



3 • - 


j-L* j 


2 * 


J = 


j l + j 2 


f L = 


E 

j| = L 


J! 


(y - 


V 


With this abbreviation, (18) takes on the form 


(14) 


I 


2 


SS F(y) e- pf(y) dy 
C sdp 


(15) 
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Xo obtain an asymptotic expansion, F(y) and f (y) are rewritten in 
;Cl s of a Taylor’s expansion about the saddle point y . 

n 


F(y) « I F 
K“1 


n 


f(y) = Z f 
K=1 


>c exponential term is exp nded so that 


( 16 ) 


*of(y) _ -pf -p/z f 9 {1 - p[f« + f y . + ...] + 


= e o e 


(17) 


n 


[f 0 + t, + ...J +...+ [£3 + + • • * ] + * • *] 


n 


:f one lets u “ n Q = % and V - V q - f)» this manipulation yields a solution 
,0 (15) in terms of a summation of integrals of the type 


00 co .2 2 

f , r m n -p /2 [f £ + 2 f £r) 4- f r| ] 

* / / £ rj e uu’ uv . VV J dt,dn 


nm 


(18) 


—00 —00 


were the limits to +°° have been taken as an approximation to . 
1 solution to integrals of the type in (18) has been found to be 


I = n + m.m 
nm 1 — 


(n + m) /2 + 1 


’fie re 


r c 


m + n 


C , = «S 

n 4- m,m 


4- 1) , tn - n . n 

£ K 2 } Z R(n,k) 

uu , * 

k = 0 


; for n + m even 


0 


; for n + m odd 
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and 


r Q 


R(n,k) 


- j 


(f f - f 

uu VV uV 


• j(n + k)/2 rC 5 ^- 


+D rc 15 -^ +i) 

(19) 


; for m + k and n - k even 
0 ; for m + k or n - k odd 


To obtain an asymptotic series , the expansion of integral (13) is 

combined in terms of increasing powers of 1/p. 

- n l “l 

If there are polynomials of £ and t| associated with the Taylor's 

n 2 m 2 

expansion of F(y) and polynomials of £ and r) associated with the 

expansion of , from the form of (19) it can be seen that these 

. - ,1, (n. + n~ + m_ + m )/2 + 1 , 

polynomials have a contribution of (— ) 12 12 . ihe 

P 

subscript J = is now associated with the Jth term of the Taylor^ 

expansion, and a value E& " n z m 2 i- s associated by Table 1 with products 

f (y) 

of f in the expansion of e 
Lt 

There is also a contribution of (p) * dire'ctly from the expansion of 
e in (17). A rule can now be formulated for generating the nth 

term of the asymptotic expansion. Specifically, one relates the nth term 
of the asymptotic expansion through terms of f^'s and F^’s for which 


,J + ZZ , 

( — ^ + 1 - k) = n 


( 20 ) 


Those terms associated with n = 1, 2, 3 are summarized in Table 2, 
Furthermore, the contribution by these combinations can be written 

-pf 


0 (-D k 

n k! 


F J f Ll • ' - f Ln C J + Zl, m + m 2 


( 21 ) 


where ~ L. + . . . + L , J + £& = m. + m_ + n. + n and by C T , , 

1 n 1 2 1 2 ; J + Ex, + 


m. 


is meant, the integration constant in (19) associated with the polynominals 
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-of (y) 

TABLE 1. ORGANIZATION OF TERMS FROM e 


I • J •«* N Corresponds to f^. fj 



3 4 56789 |l 


34 

35 

36 

43 

44 

45 


53 

54 



55 



333 


TABLE 2 ♦ VALUES OF J , k FOR 
FIRST 3 TERMS OF ASYMPTOTIC EXPANSION 




2D 


n *f n m + m 

£ D This value of C, though being different for the 

various combinations involved in F , f L> will be abbreviated to + 
The first three terms of the asymptotic expansion then becomes: 


1 “ pf n 

- e F C 
p o o 


1 ~ pf n 

e (F 0 C 0 - F n f„C, - F f C. 
2 22 134 q44 

P 


~r Ff,f.CJ 
21 o 3 3 6 


h e 0 {F 4 C 4 - F 3 f 3 C 6 - F 2 f 4 C 6 
P 


" F l f 5 C 6 " F o f 6 C 6 + 2! F 2 f 3 f 3 C 8 


+ 2! 2 F l f 4 f 3 C 8 + 21 2 F o f 3 f 5 C 8 


+ 21 F 0 f 4 f 4 C 8 ‘ 3! F l f 3 f 3 f 3 C 10 


- —■ 3 F f_f.f.C in + 77 F f 1 f XC ,} 
3! o 3 3 4 10 4! o 3 3 3 3 12 


( 22 ) 


This form of the asymptotic expansion for the first three terms 
has been programmed on a computer. Considerable effort must be spent 
obtaining the many partial derivatives needed. In fact, since this is 
a significant difficulty in using higher order terms, some thought should 
be given to developing a computer program to perform this algebra. 


4. Applications to the Free Space Problem 

The asymptotic techniques developed in the proceeding sections were 
applied to- a dipole current source in free space as a test to the approxi- 


mation. 



21 


From Maxwell’s equations 


curl E(r) = -j< 0 MH(r) 
curl H(r) = jtoeE + J(r) 


( 23 ) 


One obtains 


curl curl E(y) « <oyeE - jwyJ(r) 

E is now considered to be a plane wave and J (r) a source in the 
z direction only. In the transform domain, then 


n 2 /r 2 . 2 . 2.., ~ 1 ~ n 2 2 , 

[k -.(€ + n + ; )] e = — j(k - c) 


where 


00 + j (£x + r)Y + c z ) 

E (k)' = / jf / E e dxdydz 


and J is the Fourier transform of J(z). It follows then that 
1 


co 2 2 “jk * r 

E = ^ L. / J f d \ 

h Q 3 jwe J J J , 2 . 2 .2. 

8tt j -» k - (£ + n + £ ) 


A change of variables is now performed 


t, - r cos 3 

X] = r sin 3 sin a 
£ = F sin 3 cos a 
With this change 


x = r cos 0 
y - r sin 0 sin <J> 
z = r sin 0 cos <j> 


(24) 


~ , 2 ,-,2 2^ „ -jTr cos y 

F sin 3 e dT da d3 (25) 

C(a)c(B) k^ - r‘ 


E = —■ ~r— f f f J , , 

oTT J 0J6- j-. / \ * / a n i ^ 


where cos y = cos 0 cos 3 + sin 0 sin 3 (cos(a - $)) and C(a) , C(3) are 
paths in the complex a and 3 plane. 

The f integration is done by a method of contour integration resulting 


is 


E(z) = 


~ ~ -jkr cos y 

/ / J sin 3 e dad0 


(26) 


toe tt c(a) C(3) 
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An asymptotic expansion of the integral in equation (26) is now used. 

We have, using previous notation, 

F(ct, 3) = J sin^ $ 

f(a, 3) = cos y (27) 

and p - -jkr 


And the saddle point is at a = <p , 3 = 0« 

Let us consider a sinusoidal current distribution 
J = (x,y,z) = 1^ sin (k(H - |z|) f° r which the field is already known. 

This is the form for the current distribution on a thin dipole of length 

{43 


2H. The solution' 


E(z) = ~j 301 (- 


for E is 
z 


-j3k n 


-j0R 


+ 


R, 


2 2cos(3H) 

n ^ / 


(28) 


'1 2 

The geometry of this dipole is shown in Figure 1. 

For this current an asymptotic approximation to (26) can be found, 
and the results compared with (28). Moreover the exact solution shown 
(28) can be expanded in terms of a series of (-;) and compared with the 
asymptotic expansion, which is also written in terms of a ^ series. 

It should be noted that the calculation of the first three terms 
s ~of the asymptotic expansion is facilitated by the condition f(<£, 0) = 0. 


2H 


1 X_ 
9* 10 


The approximate and general solutions were generated for — r- = — , 

A o 

and the radius of the dipole was taken much smaller than a wave- 


length. The first three terms of the asymptotic approximation agree 
exactly with the first three terms of the expansion of the exact solution 
(28) in terms of a series of (“) . At a distance of .4 wavelengths, the 
asymptotic solution was within 1% of the exact solution, and at .21 the 
asymptotic solution was within 10%. 
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Within . 2X the value for the approximate solution quickly diverges 
from the exact solution. Thus, it appears, and this was expected, that 
for very small values of 2 , E(z) cannot be evaluted using an asymptotic 
approximation. However, the asymptotic approximation using the first 
three terms is quite good for intermediate and far distances* 


5. Dispersion Surfaces 


As mentioned earlier, to find the approximation to the field at an 
observation point r one expands about the saddle points, which correspond 
to those points for which dk • r = 0, In physical terms, one is finding 
the k vector of plane waves which contribute mostly to the far-field in 
a particular direction. These correspond to those points k for which 
the perpendicular to the dispersion relation is in the r direction. 

The geometric problem associated with saddle point integration can 
be stated, therefore, as given an observation direction r, generate the 
k vectors for dk • r = 0 and the dispersion relation being satisfied. 

The dispersion relation can be described 1 using Maxwell’s equations 
for a plane wave 



(29) 


So Maxwell's equation becomes 


kxE — wB 
kxH - -o)D 


Letting K=wyee=K one has 
00 


kxkxE - -KE 

(k k T - k 2 + K) E = 0 


(30) 


(31) 

(32) 
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For non-trival solutions, therefore, 


det (kk T - k 2 + K.) = 0 


(33) 


which describes the dispersion relation for the media. This dispersion 
relation behaves differently for different types of media. 

For an isotropic space (33) implies that 


2 2 
k = w \i e 
o 


From which it follows that 


2k dk = 0 


(34) 


(35) 


One conclusion of (35) is that in isotropic regions k and the observation 

angle are in the same direction at the saddle point. 

Consider now the uniaxial plasma. This type of region corresponds 
to a plasma with a large external magnetic field with an E matrix of 


the type 


e = 


e A 0 0 

0 0 

0 0 e. 


(36) 


One solution to the dispersion relation is that 


2 2 


And for this value of k 


(37) 


2k dk = 0 


( 38 ) 


so that k again is in the same direction as the observation angle at 
the saddle point. However, other solutions to the dispersion relation 


are available. 
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Letting P = (k 2 - k) X kk^ (33) becomes 


det (I - P) E = 0 (39) 

Because the image of the operator P is always d, P has rank one. Thus, 


Or 


det (I - P) = 1 - trace P = 0 


trace P = 1 

T 2 -1 

k (k - K) k = 1 


(40) 


(41) 


So, recalling that dk * r = 0 at the saddle point, (41) can be 
utilized to obtain the r direction by taking the differential with 
respect to k. Sepcifically , 


T 2 -1 T 2 -1 

0 = dk (k - K) k + k (k - K) dk 


T 2 _ i 't o -1 

-2 k (k - K) 2k k" (k - K) 


(42) 


0 = dk T [(k 2 - K)' 1 k + (k 2 - K T )“ 1 k -2 k T (k 2 - K)^ 1 (k 2 - K)" 1 k kj (43) 


Thus, the observation point is in the direction of a vector k where 

R = (k 2 - K)" 1 k + (k 2 - kV 1 k -2k T (k 2 - K)~ X (k 2 - K)" 1 k k (44) 

Equation (44) holds for the general lossy anisotropic plasma also. 

It is necessary that an inverse relationship be found. That is, 
given a radial direction r, a direction of k must be found, k may then 
be scaled to fit the dispersion relation (33). 

6. Uniaxial Case 


The problem of determining the direction of k at the saddle point 

T 

has been solved exactly for the case K = K . Equation (44) reduces to 

1 „ ,,2 t , n -1 , , T /. 2 „.-lT .. 2 „v~l , , 

— R - (k - K) k - k (k - K) tk - K) kk 


( 45 ) 
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Some physical interpretation can be given to (45). Recall that 

T 2 2 

(k k - k + K) E •“ 0 so that for (k - K) nonsingular 


E = (k 2 - K) 1 k k T E (46) 

° r 2 -1 

E = (k * E) (k - K) X k (47) 

12 11 
Thus ~R - CE - 2C (E • E) k for C - ~ - ~ , and -^R = CE x (k x CE) 

2 

= C m (ExH) (48) 


Equation (48) has an interesting physical significance. It implies that 
the normal to the dispersion surface is not in the direction of the 
Poynting vector, the direction of energy flow. 

Now some dual relations can be constructed to solve the geometric 
problem. Let 

s= or^kT R . . (49 > 

Then 

S X B = - SX (k X E) 

0) 

S X B = - £ (S * E) k - (S * k) E] (50) 

S X B = -E 

Similarly 

S X D = H (51) 

A dispersion relation for S may now be derived. 


SX (S X D) I ~ 1 D 

\i£ 


(52) 


1 Z -1 T 

Letting V = £ , V = V for the uniaxial case. So, 

PE o 


(SS T - + V) 


D = 0 


(53) 
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This implies that for non-trivial solution 

det (SS T - S 2 + V) = 0 (54) 

so that 

det (1 - (S 2 - V)" 1 SS T ) = 0 (55) 

and 

S T (s 2 - v)" 1 S = 1 (56) 

Taking the differential of (56) with respect to S results in 

[dS T ( (S 2 r V) _1 S - S T (S 2 - V)' 1 (S 2 - V)" 1 S S)) = 0 (57) 

The results of (57) must now be examined. From (53) we have that 

D = (S * D) (S 2 - V)” 1 S (58) 

From which it can be shown that k f — (S — V) S — S (S — V) (S — V) S S 
1 

is in the D X H or k direction. 

Thus for a medium in which z = Z one is able to obtain the direction 
of k at the saddle point given the direction of the observation point. 

This value of k is then scaled to fit the dispersion relation (41) . 


7. General Magnetoplasma Region 


For the general anisotropic region, £ has the form 


£ - 


"'i 


“je 


I 0 , 0 I 

One has not been able to find an analytical solution to the geometric 
problem of obtaining the k vector for which dk • r = 0. 


( 59 ) 
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Optimization techniques have been tried in the anisotropic region. 
However, the general lossy anisotropic region will have many multiple 
solutions for which a searching scheme must be devised. 

8. Conclusions 

Using asymptotic techniques an approximation to the solution of the 
ik * £ 

integral / F(k) e dk can be made for F(k) varying "slowly" at a 

C 

reasonable magnitude of (k • r) . This approximation, expands this integral 
about a saddle point. This has been done for the one— and two-dimensional 
problem. Preliminary results indicate that in free space, the estimate 
is with 10% of the exact solution at r = .2 wavelengths and within 1% 
beyond r = .4 wavelengths. 

The problem of finding the saddle points corresponds to finding the 
k vectors which contribute to the field in a given direction. An analytic 
solution to this problem has been presented for the special case in which 
the dielectric matrix has the symmetry e = e T . Otherwise, optimization 
techniques are suggested. 
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II. 3. Laboratory Measurements 

a. Impedance Measurements - Effect of Power Level and Voltage Bias 

In measuring the variations of antenna impedance with the input power 
level some difficulties with the vacuum system were observed. The following 
modifications improved the consistency of the results considerably. First 
the vacuum system was thoroughly cleaned and new gaskets were installed. 

The diffusion pump jet assembly was found damaged due to oxidation of the 
diffusion pump oil. A new jet was then installed and silicon oil was used 
to prevent further accidental oxidation. The helium flask was replaced 
also. A second mechanical pump was added to the original system and 
connected to the backfill port of the pumping station. This pump provides 
the capability to continuously exhaust the plasma tube without the fear of 
damaging the diffusion pump by accidental breakage of vacuum or other 
contaminations. It also keeps the plasma tube at low pressure even when 
the tube is not in use. This reduces the pump down time considerably. 

An electronic discharge current regulator was built to stablize the current 
and it worked beautifully. The RF system was revised also. When operating 
at 300 MHz a RF amplifier was added to boost the available power to about 
2 watts. With the addition of power attenuators and a directional coupler 
to a power meter it is possible to vary the power level from 2 watts to 
less than 5 milliwatts and to continuously monitor the input power level. 

At 800 MHz and above, the oscillator/amplifier combination is replaced 
by a powerful oscillator at our disposal, while the remainder of the RF 
system remains the same. Due to the high power level involved in the 
measurements it was found necessary to place most RF equipments inside 
an anechoic chamber in order to avoid erratic readings on VSWR meters. 
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Many efforts have been expended to determine the impedance variation 
of the antenna when the input power level was varied. Two frequencies 
were used, 800 MHz and 1 GHz. Several typical results are given in 
Figures 1 to 3. It is interesting to note that all the impedance varia- 
tions as the power level increases follow the same general pattern. This 
behavior is not completely understood at the present moment, but the 
results are very consistent. The dipole started to glow when the input 
power level reached few hundred milliwatts. The glowing is an indication 
of induced ionization by strong RF electric field near the dipole. We 
expect that the electron density increases under this circumstance. This 
is probably the reason why the impedance curve bends backwards near this 
region as it does when the electron density is raised by increasing the 
discharge current. At low power we noticed that even for a level of a 
few milliwatts the impedance may change drastically. In Figure 1, points 
1 and 2 correspond to power level of 1.6 mW and 7.3 mW respectively , but 
there is a substantial difference in the impedance. Similarly in Figure 2 
points 4 and 5 correspond to the power level of 5.0 mW and 10 mW respectively. 
However, the variation of impedance from point 1 to point 3 may be con- 
sidered as small. Their power levels are 79 pW and 1.3 pW respectively. 
Therefore from these measurements it seems that 1 mW could be considered 
as the threshold for small signal approximation. Since in practice 
measurements are commonly performed at a power level of several milliwatts 
or higher, this casts a great doubt in the validity of many existing 
results. Therefore it is concluded that any antenna impedance measurement 
without reference to the input power level may be meaningless. Another 
interesting feature is that when the external magnetic field is zero, the 
impedance of the antenna shows little changes as seen in Figure 3. 


r 
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b. Observation, of Nonlinear Effects 
Harmonics Generation 

The generation of second harmonic in plasma has been observed with 

a spectrum analyzer to monitor the reflected power from the antenna and 

also by tuning a receiver to the second harmonic frequency. The second 

harmonic power level was measured against the magnetizing field for 

various bias voltages. The results are shown in Figure 4. In all cases, 

a local maximum is shown at the cyclotron resonance. When the bias 

voltages were varied, the more positive the bias is, the higher is the 

second harmonic level as might be expecced. However, the second harmonic 

level eventually became saturated, when the bias reached several tens of , 

volts. Under certain circumstances, the third harmonic has also been 

observed but no attempt was made to measure the third harmonic level. 

When modulated signals were fed into the antenna, the modulations were 
* 

detected by a receiver at the second harmonic. 

Luxemburg Effect 

In this experiment two transmitting antennas in the plasma were used. 

One was fed with a modulated signal with carrier frequency f^ at a fairly 
high power level and the other with a continuous wave at frequency f^, A 
receiving antenna was located several feet away from the plasma tube and 
the receiver was tuned to f^* It was found that the modulation or audio 
signal was detected through a speaker, indicating that part of the modulation 
was transferred to the wave at When one of the two transmitting antennas 

was placed outside of the plasma chamber. A similar phenomenon was observed 
but with weaker modulation. This is the so called Luxemburg effect or 
crossraodulation effect. In ionosphere it is believed that the electron 
heating is the main cause of crossinodulation. In a laboratory plasma many 
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other nonlinear effects could produce crossmodulation between signals. 
It is not clear which nonlinear effect is the dominant contributor. 
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II. 4. Reflection and Transmission of Waves in an Electron-Heated 
Nonlinear Plasma 


X . Intro duction 

There are various kinds of nonlinear phenomena in a plasma. Tsytovich 
(1) has made an extensive and detailed discussion on the nonlinearity due 
to the (v • V) v and v x B terms in fluid equations whereas Whitmer and 
Barrett (2-5) considered the same set of equations and calculated the re- 
flected and transmitted v7aves , including harmonics, through a plasma 
layer. Keller and Hillman (6) assumed a very general constitutive re- 
lation with nc particular reference to physical processes and solved the 
Maxwell’s equations for several propagation problems by using perturbation 
method. Ginzburg (7,3) discussed the nonlinearity due to electron heating, 
in particular its application to the cross modulation of waves in ionosphere, 
the so-called Luxemburg effect. In his elementary theory the state of 

the plasma was characterized by two variables, the average electron velocity 

■+ 1 ^ 
v and the effective electron tecoerature T . Ginzburg shewed that ir the 

e 

applied field intensity is smaller than the "plasma field" the state 

of the plasma changes slightly under the influence of the field . however 

when the applied field intensity becomes larger then E , the parameters 

P 

of the plasma change significantly. Since in commonly encountered plasmas 

the value of E is small, significant nonlinear effects may result even 
P 

for moderately large applied fields. In this paper, the nonlinear effect on 

the reflection coefficient for a plane wave incident on air-plasma interface is 
investigated . The reflection from a semi-infinite nonlinear plasma is 
considered first and from a conducting plate backed plasma layer next. 

The reflection coefficient as a function of incident field strength is 


solved for both cases. 
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II, General lo rm.ul a tion 

In a nonlinear plasma, the mediur; property depends on the applied 

field strength. Ginzburg (7) has shown that if the nonlinearity of the- 

plasma is mainly due to electron heating processes, the local effective 

temperature of the plasma is a function of the magnitude of the field 

strength E . At low frecueiicies the properties of the plasma are modulated 
v 

by the applied field and the problem becomes rather complicated. Therefore 
we will confine cur discussion to high frequency limit such that of » v. 

In the high frequency limit, the field varies much faster than energy re- 
laxation time , therefore we can replace E by its rms value E or aside from 
an unimportant constant, by its maximum value E^. from Ginzburg’s elemen- 
tary theory the relation between T and E can be expressed as follows 
J J e o 


2 2 

0 T E n 0) + V 

T 2 H =|*(/) 2 - J — 5 - e — = 1 + (^) 2 if v 1 

p CO + V (T^) P 


(I) 


where T is the equilibrium electron temperature of the plasma in tne absence 

of external field, v and v are effective collision frequencies in the 

o 

absence and in the presence of external field respectively. V can be 

considered as a function of T or E . E is the so called "characteristic 

e o p 


plasma field ; 


E = /3kT 
P 


hr 5(a* 4- V “) 
1 o 


( 2 ) 


where 6 is the average relative fraction of energy transfered during a 

collision between electrons and heavy particles. It is seen that when 

E << E the heating effect will be small and T - T, but when E > E the 
op e o — p 

electron heating processes becomes significant. It is convenient to 
introduce the concept of dielectric permittivity and conductivity for a 
plasma. It is well known that 
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LO VW X 

£ = 1 - . E c « p- (3) 

2 1 2’ // 2,2. • v ' 
a) + v 4tt(w + v ) 

where is the plasma frequency. Since v is a function of E q , there- 
fore e and o are also functions of E 

o 

Consider a one dimensional case, the electric field inside the plasma 
satisfies the simple wave equation 


d E , 2 , E 

■ — r + 0) £' — r = 0 
dz c 

and here e f = E - i = (n - iX) ^ 


- n -i i ‘ 2 ,4710.2 . , 

In plasma, usually c » ) is satisfied 
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W 

(5) 


(6) 


We will also assume that collisions are mainly between molecules and 
electrons. The collision cross-section is approximately a constant in- 
dependent of veloicty such that V = TV^. In high frequency limit X “ TX c 

« 2 
where 2ire ftv 

X - ■p^rrr-: 

° I 2 ■ 

M- <hs_ 

and as a first- approximation n is a constant. In this case we can use the 
geometrical optics approximation. A solution for the wave equation (4) 
can be written as 


E *= A exp[io)t + i “ / dz]- E exp [ i (cot - — f Z ndz)] (7) 

C o c 0 

For simplicity we denote as the magnitude of the field. Sometimes it 
is easier to use a differential equation for instead of (7). Differentiation 

of E with respect to z results in 

o 
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d.E 

f: - - - X(E ) E 
dz c o o 


( 8 ) 


Now equations (1) and (8) are to be solved simultaneously subject to 
suitable boundary conditions. 


III. Reflection from a Semi-Infinite Plasma 


First consider a simple problem, a plane wave incident normally cn 
a semi-infinite plasma which is assumed to be uniform and isotropic. Let 
the interface between free space and plasma be tha z ~ o plane, and the 
electric field is in the x direction as shown in Figure 1. 

Substituting >' into equation (8) and make use of equation (1) 

- K 

„ 2 


dE 


o 

dz 


to 

— TX 

C A 0 0 


E — — — ^ — — y (t — 1) T 
& c o 
o 


( 3 ) 


dE 
o 

dz 


z=0 


E 2 ' 

“ - EW c X o T o (T o - ’ T o = T(Z=0) 

o 


( 10 ) 


From Maxwell ’s equations the magnetic field can ce computed 


iU r, / n s 8E 

B (z = 0) = u 

c y dz 


dE 

0 

i (wt - — f Z ndz) 

p A* A 

dz 

C L U 

z-0 

z=0 


z-0 


, . w d ,z , . x(tut / ndz) 

+E (z=C) [-i — — — / ndz J e c o 

o c d z o * 

z~0 


(ID 



J-/ 4.. 


.. >2 

y ' o"" “ ■ ' toe v ] 


B (0) = /£ E (0) [1 + (• 


5.(0 1 

e Y e 




( 12 ) 


*= C E o (0) e 1 ^ c iut 


where £ = 


. -1 ,2iro 

“ tan ( t ) 

we o 


(13) 


The reflection coefficient R and transmission coefficient T then can be 
calculated easily by Batching the boundary conditions at z = 0. 

It can be shown that 


j E | ‘= \Lz Ce 1 * 
h + ce 1 * 


and as usual 1 + R = T. 

Figure. 2 shows the reflection coefficient and transmission 
coefficient as functions of incident field intensity. The horizontal 
axis are incident field intensity normalized with respect to E . me 

parameters chosen were typical for a laboratory plasma and the^harac*. eristic 
plasma fields E p are of the order of 100 volts/cm. 

^^S-jsS£££l£^_f£ca_a_Jioaline2r^ Pla S3a Laver Bac ; vrf , p .^ 

■ Now consider a plane wave incident normally onto a nonlinear plasma 

slab with thickness d and backed by a perfect conducting plate as shewn 
In Figure 3. 

It has been shown that when a? » v* n can be considered as a constant. 

In the present case the incident, field Ej _ and the reflected field E „ i r . 
the plasna can then be expressed in the following form 

r 

1 


El = ae 


-ikz - - / z piJ z 

£ C O 


V E 2 = be 


^ p ~ S Z PA 2 

e c- n * 


( 14 ) 


where k - and for simplicity the time factor e 1 ** has bee, suppressed 
One of the difficulties of solving a nonlinear problem is the in applicable 
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of the superposition principle. For the pits eat problem the absqrp oa 
coefficient X depends on the total electric field strengths which is the 
sum of lh and f ' ; therefore the unknown functions E , E and X have t0 

12 i. ^ 

be solved simultaneously. The equations are 


r 


d ! E ii 


+ 7X (H 0 ) \\\ = 0 


■N 


} dll] 

; 2 . - £ v-(e ) ieJ 

] dz c X ° • 2 1 


= 0 


(15) 


(16) 


t 2 , 

nr i = 


(17) 


where E = E 
o o 


From (14) and the boundary condition for E^st z = d, the biaguitude of 
can be written as 


*• ; ■» 


E = 1 E + E ! = y ( 
o « L 1 c 2 [ 


4 — / 2 Yd~ 2 — fZ Ydz , xt ‘ 

(z) [1 + s c J d X ^ - 2 cos 2k (s - d),e. c ' c X • 


(13 


where E (z) = j (n) ] * X is a function of E o , 


Introduce a new function U(z) = 2 — x /j then substitution of 

C O c 


equation (IS) into (17) yields 

2 2 

2 r E n ^ " i} 

~ 2 -b _ 2 

E a ' . 2U . U 


(19) 


[1 + e - 2e cos 2k (z-d) ] 


du to 

where E . - E (z=d). Making use of ----- - 2 — x T, equation (19) 
d u/ c o 

converted into a nonlinear differential equation 


can ce 


,22 

4uj x 

— j » [cosh U (z) - cos 2k(z-d) ) -1- 


„ , n 2 2 2 
Jtt \ 2 tZ. a) x 
dU ! d ,v o r 


( 20 ) 
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It appears that this nonlinear differential equation cannot be solved 

in closed form, and therefore numerical techniques are used to determine U and 
dU 

dz* ^ ne set calculated results are given in Figure 4 for E^ nc /K ~ 0.18, 
0.47, 1.16, 3.29, 13.0 in which the total field strength E q is plotted 
against z. The plasma parameters chosen in these calculations are the same 
as before and the slab is five wavelengths thick, ’Then the fields is weak 
as in 4(A) the solution is essentially the same as that for a linear plasma 
slab; i.e., the heating effect is insignificant. However as the field in- 
tensity increases the absorption also increases and the field distribution 
inside the plasma alsb changes considerably, becoming strongest near the 
interface. In contrast, for a linear plasma the shape of the field dis- 
tribution remain unaltered, regardless of incident field strength. The 
reflection coefficient |r| is plotted against the normalized field intensity 
in Figure 5. It is interesting to note the existence of a minimum 

reflection fcr a particular value of E. 

me 

V . C o n c lu si on 


'The reflection coefficients Jk{ from a semi-infinite plasma and a 
plasma layer backed by a conducting plane are determined as a function of 
incident field strength. In the first case it is found that | R J is in- 
creasing mono Conically and approaches 1 as the incident field increases 
touard infinity (within the range of validity of our theory) . For the 
second case, the reflection coefficient has been found to have a minimum 
for a certain value of incident field. For large value of incident field 
the reflection coefficient increases steadily and approaches 1 as in case 
one. The total field strength inside the plasma layer is evaluated also. 

It is expected that tnis particular nonlinear effect can be observed in the 
laboratory by measuring the reflection coefficient from a discharge. • For 
ionosphere the effect may be too s frail to be significant. 
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Figure 1, Wave reflection from a serai -in finite plasma. 

Figure 2. Reflection and transmission coefficients as functions of 

incident field intensity. 

in. q 

f - 6 x 10 Hz V = 1 x 10 /sec 

o 

N - 1 x 10 12 /cm 3 6 = 1 x 10 -3 

T = 1 x 10 4 °K 

Figure 3. .Reflection from a plasma layer backed by conducting plane. 


Figure 4 (A) -4(E). The variation of field intensities inside the plasma 


■layer as the incident field increases. At the air-plasma 


interface E. / E equals: 
inc p 


' (A) 

E. /E = 0.18 
me p 

<B) 

E. /E = 0.47 
inc p 

(C) 

E. / E = 1.16 

anc' p 

CD) 

E. /E » 3.29 
inc p 

00 

E. /E = 13.0 
me p 


Figure 5. 


Reflection coefficients as functions of E. / E . 

me p 
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• 5* Profile Inversion of Inhomogeneous Media 

1 . Background and General Description Of the Proh 1 era 

We have investigated the problem of determining the electrical permittivity 
of a medium which is uniform in the x and v directions and varies ac- 
cording to the profile function, ’<(/.), in the z direction. The mathe- 
matical problem associated with this profile inversion problem is typical 
of many other remote sensing experiments. 

1 2 

Using the techniques of Gelfand and Levitan and of Harchenko, many 

3 

authors, e.g., Becker and Sharps'^ have reported the solution of this 

type of problem. All of their work, however, requires that; the response 

of the medium be known as a function of the temporal frequency. It is 

also required that the medium itself be independent of frequency, i.e., 

nondispersive. Often, however, it is desirable to make measurements at 

a single frequency. ^Therefore, in this paper we consider the remote 
1 

probing problem in the angular spectrum domain. Here the independent 
variable, 0, is the spatial frequency which can be associated with k sin 0 
where k is the free space wave number and 0 is the angle of incidence of 
a probing plane wave. 

The technique by which the response of the medium is measured is to 
illuminate the medium using a spatially confined source. The total field 
.at the interface is then measured. The information is then converted to 
the angular spectrum domain by taking the Fourier transform with respect 
to the x direction. 

In the next section, the formulation of the problem is presented. The 
remaining sections are devoted to techniques for solving the problem. 
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2. Formulating the Problem . 

The geometry of the problem to be considered is shown in Figure 1, 
Note that there are nc variations with respect to the y direction. The 
dielectric is also uniform in the x direction and varies in the z direction, 
as k(z)e^ where is the permittivity of free space. The permeability 
is assumed to be constant and equal to that of free space. It is assumed 
that a known electromagnetic wave of a single frequency is incident from 
the left on the air-dielectric interface at z = 0 and that the tangential 
fields at the surface of the dielectric can be measured for all values of 
x. Of course, the actual measurement could be performed at some plane 
z = a, ot < 0 , and then the fields at the interface could be computed since 
the propagation in free space is completely known. 

If the incident zield is polarized such that the components E , E , 

X z 

and are all zero, then propagation in the medium is completely described 
by the scalar equation 


2 2 

V v(x,z) + to p K(z)e v (x , z ) = 0 
o o 


( 1 ) 


(to is the angular frequency of the incident wave). The function v(x,z) 
may be identified with E or H , which ever is more convenient for the 

y 2 

problem at hand. 

To obtain the differential equation in the angular spectrum domain, 

we introduce the Fourier transform relations: 

00 

u($,s) = / v(x,z) exp(-i0x) dx (2a) 


v(x,z) = 2^ in u (3»z) exp(i0x) d0 


(2b) 


Then Equation (1) becomes 


+ [u 2 y o) c(z)e o - B 2 ] u(g, z ) = 0 
dz 


( 3 ) 



;>4 

This is one form of the differential equation in the angular spectrum 
domain . Recall that to, p , and £ q are constants, 3 is the independent 
(spectral) variable, and k(z) is the unknown profile function. 

. For the work which follows, it is convenient to introduce a different, 
but equivalent, form of (3). That is. 


where 


+ y 2 u (Y,z) = Q(z)u(Y,z) 

dz 


2 i 2 a 2 
Y = k - p 


Q(z) - -k [*c(z) - 13 


2 2 
k = W P £ 

O D 


( 4 ) 


Note that ‘if we determine Q(z), k(z) can easily be found. Furthermore, 
note that we intend to apply (4) in the source-free region 0 <_ z < 05 . 

.3,, The Solutions <fo(Y>z) and f(Y> 2 ) 


■Equation (4) turns out to be of the same form as the one-dimensional 
Schroedlnger wave equation. This equation received considerable attention 

in the 1950 f s from physicists who were interested in scattering by 

potentials. Host of the work on the inverse scattering problem has been 

1 2 

based on the techniques of Gelf and-Levitau and Marchenko , An extensive 

review of work on the inverse scattering problem has been presented by 

_ i , 4 

Faddeyev. 

We will assume that the profile function Q(z) is real and bounded 
and that 

/ z I Q (z) [ dz < 00 (5) 

o 


r 
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We then consider two independent solutions of (4), $( y,z) and f(y,z), 
satisfying the boundary conditions. 


<Ky,0) = 0 


d<j>(V,0) . , 

dz 


(6a) 

(6b) 


and 


lim 

Z-KO 


[f (Y,z) 


-iyz 

e 


] 


1 


(7) 


The solution f(y,z) corresponds to the situation shown in Figure 1. The 
solution 4>(y,z) corresponds to situation where a perfectly conducting wall 
is placed at z = 0 and the incident wave comes from z = + It is de- 
monstrated in Faddeyev's paper that 


2 ^ <J>(Y»z) = sin[y'z - ^(y)3 


( 8 ) 


i/he re 


f(Y.O) - A(y) e 


i^(y) 


(9) 


4. The Inversion Technique 

. 5 

Mittra, et al . , have introduced the following representation 
of f(y,z). 

f<Y.«) = - P (i YZ ) + ^ C «n.*) 

In this case, the profile function Q(z) is given by 


( 10 ) 


Q(z 


) = ~ /°° [i - n f (n»z)] exp(irjz)p (n) d 

TTl — 00 dz 


on 


( 11 ) 


The function P(n) can be determined from f(y,o) by considering (10) at 

2 = 0 , 

f ( y>o) = 1 + 2 ?I £> f (n,°) an (12) 
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These three equations provide a technique for solving the. inverse scattering 
problem in the spectral domain. That is, a technique for obtaining Q(z) 
from a knowledge of f(y,o). 

It is common for inverse problems in mathematical physics to be ill- 
conditioned, This means that small errors in the data result in large 
errors in the computed profile function. Mathematically, ill-conditioning 
is usually associated with the solution of a first kind Fredholm equation. 

We note that (12) is indeed a first kind equation for the unknown function 
P(ri) , Thus we v/ould expect to have trouble solving (12). On the other 
hand, once P(n) is specified, (10) becomes a second kind equation for 
f(Y» z )~' 2 is fixed and y is the independent variable. We, therefore, 

would not expect (10) to be ill-conditioned. Thus, we are led to the 
conclusion that if we can obtain P(r|) accurately from f(y,o), we will have 
a procedure for obtaining Q(z) which is numerically stable. 

We have been able to show an analogy between (10) and che representation 
which Marchenko uses for f(y,z). This has allowed us to determine that 
P(ri) is given by 


P(n) = 1 - e' 


-imn; 


Using (13), we then solve (10) for f(y,z.) for several values z » The 

J j 

profile function Q(z) is then obtained from (11). 

It should be noted that (13) is valid only if f(y,0) has no poles 
or zeros in the upper half plane. This condition implies the absence of 
surface wave solutions to (4) . If surface waves are present in the 
problem, the inversion must be carried out in two steps. The procedure 
is exactly analogous to the method described in Faddeyev's paper. 


r 
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One further point to note about this inversion technique is that P(rj) 
as given by (13) utilizes only the phase function $(r|) . The condition under 
which (13) is valid imply that f(o,rj) is an analytic function of rj in the 
upper half plane. Thus a knowledge of (n) is sufficient to determine A (n) , 
and vice-versa. Hence, it is sufficient to specify either the magnitude or 
the phase for this case. 

. • 5 . Numeri cal Solution of Equation (10) and Example s 


Using the technique described in the previous section, several 
examples of profile inversion have been solved. To do this, (10) must 
be converted to a matrix equation which can be handled numerically. 


Before proceeding, with the numerical formulation it is convenient to 
multiply (10) by P(y) and then, defining \b( y,z) £ P(y)- f(y,z),.we have 

*<Y.z> - IS 1 L *<n.s) dn = p(Y)e iyz . (14) 


Investigations of the asymptotic behavior of f(y,z) and P(y) indicate 
that as Iy| -> CO 

P(Y) ^ 0(i) (15) 

f(Y,z) ^ 0(e 1YZ ) (16) 

Also, as z «>, 

f(y,z) ^ 0 (e 1 ^ 2 ) ( 17 ) 

Furthermore, if the medium consists of several homogeneous layers, then 
f(YjO) will not have any branch singularities. It can also be shown that 


f(-y,z) = f (y»z) (18) 

where the bar denotes complex conjugate. These criteria have led to an 
approximation of ^(y,z) of the form 


r 
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ifi(Y.z) = 


i.yz v 

e 1 £ 


11 * B j(z)v 2 + i C , (a)y + ( 2 ) 


j=l 


4 4 

Y + P. 
J 


(19) 


(This form is not the most general rational expression that could be used, 
but the results obtained using this form have been quite acceptable.) 

We are now ready to convert (14) to a matrix equation. First, we 
take the inner product of (14) with the functions {6(y - This 

generates m equations of the form 

P(Y,) „ exp[i(n + Y,>] (20) 

U\,z) - 2iT~ io 'Hi- 2 ) j f+Y dn = P(Y k )e ' k = l,2,...,m 

k 


(Recall that z is a fixed parameter in (20)). Next, we substitute (19) 
into (20) to obtain 

1A j (z)[H y - + - i li> 


+ 1 C j (z)tH kj - ^ + D j^ Z )^ H kj - ^ 


( 21 ) 


P(Yj.)e 


4 V 


where 


k = 1 , 2 , . . . , m 


^ ■ _ iY k z 




- e 


Y 


4 4 

Y + P. 

J 


£ = 0 , 1 , 2 , 3 


.£ _ P ^ Y k^ 00 n A exp[i(2n + yjz] 


/ { Li _} {. 

kj 2ni » r) + y. 


-7 r } dn 1 = 0,1. 2, 3 

0 + p. 


Equation (21) can be cast into the form of a matrix equation which is 
easily solved by numerical methods. 


■ Several comments should be made concerning this formulation of the 

problem. First, (13) is used to obtain P(y) so that the required data 

are the values of the phase function ijj(y) at the points y^j lc - 1, 2,,.., m 

l 

Second, the integrals defining the coefficients I can be evaluated by 

kj 

closing the contour in the upper half plane and using Cauchy T s Residue 
Theorem. (The pole at q = - y^ is outside the contour and does not con- 
tribute to this evaluation). The ability to represent ty(y,z) by functions 
for which the integrals can be analytically evaluated is important to the 
success of the numerical solution of this problem. Most of the examples 
we have solved required a matrix of order 16. The computer time needed to 
generate the matrix elements and to solve the matrix equation is about one 
second on a machine such as the IBM 360. 

Another point to note about this formulation, is that the substitution 
of (19) into (11) yields a summation of terms similar to the . These 
can again be evaluated using the Residue Theorem, Finally, note that the 
choice of the parameters p used in the representation of ^(Y» 2 ) is fairly 
easy. We simply note that i^(y,0) is known from the data about f(y,0). 

4i(y,0) *= P(y)f(y,0) * 

= A(Y)e ^(Y)„ A(Y)e -^(Y) ,,(22) 

= i2Im{x (y,0) }. 

Thus, we solve the equation (21) corresponding to z = 0 for a few choices 
of the p.'s and choose the set which gives the best fit to ^(y,0) . 

Figure 2 indicates the type of results which have been obtained. 

All of these examples correspond to cases where no surface wave exists. 

Note that several of the inverted profiles show large deviations near 
z « 0. This difficulty seems to be linked to the numerical behavior 
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of the equation (21). The location and the number of poles used in the 

representation (19) seems to effect the behavior of the profile more near 

i’ 0 chan for larger z. The best results are obtained when one uses the 

Eilnicaum number of poles which still gives a good approximation to the 

specific surface fields. Figures 2c, 2e, and 2f show results obtained 

v;hen the data was polluted by noise. The noise did not greatly deteriorate 

the results of Figures 2c and 2f, but the inversion of Figure 2c shows 

large deviations in the region z < 0.5. This result is not surprising 

because the profiles in Figures 2d and 2c were somewhat difficult to 

obtain even when no noise was present. This difficulty appears to be 

related to the fact that the true profile "oscillates . 11 Mittra and 
5 

Mostafavi also found that surface field data could not be used to dis- 
tinguish the true profile from an erroreous profile when the true profile 
has oscillations. 

6. Another Inversion Technique - 

The concept of transformation operators which was introduced by 
Friedrichs^ has been used as a basis for developing inversion schemes 
(See, for example, Faddeyev). The idea behind this process is that under 
appropriate conditions an operator U exists which transforms a given 
operator into another operator In our case, and are differ- 

ential operators and the transformation operator U can be determined from 
f(Y,0) in such a way that 

L 2 = U 4 'u" 1, (23) 

° operator corresponding to the profile function Q(z). In other 
•' vrf Is, if ^ and 4^ are eigenfunctions of and L 0 , respectively, then 

$2 = ^ (24) 


f 
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Considering the technique introduced in section II, we might say 
that (10) relates the eigenfunctions, f(y,z) to the free space eigenfunctions, 
exp(iyz). Using techniques similar to those used to derive (10), we have 
developed the following relationship between ^(y,z) and ^(yjz) which 
are solutions of (4) corresponding to profile functions Q^(z) anc ^ ^2^^ 
and satisfying conditions (6) . 

<f> 2 (y,z) = ^(y^z) + f° o <b 2 (r{ y z) £w x (n) - w 2 (n)J^(z,y,n)do (25) 

where 

<Kz>Y>n) = ♦ 1 (n»t)^> 1 <Y,t)dt (26) 

W.(n) - A.(n)' 2 = r 1 = 1,2. (27) 

1 1 | f i (r, »o) 1 2 

Recall that we assume the operator to be known. This means that Q^(z), 

$ (y,z) and W n (rj) are all known. Also, W_(n) is known from the given value 
of f 2 (n,o). Thus, for fixed z, (25) is a second kind Fredholm equation for 
(y^z). The equation (27) defining W, (t|) was derived by establishing an 
analogy between (25) and the Gelf and-Levitan formulation of the inverse 
problem. As was the case with the method presented earlier, we could choose 
some fixed z, say z q , where <f) 2 (y,z o ) is known and not identically zero. 

Then (25) becomes a first kind Fredholm equation for W 2 (Tl). As stated in 
section 4, regarding (12) this procedure would be less desirable than 
making use of (27). Again however, (27) is valid only when no surface wave 
solutions satisfying (6) exist. If such surface waves are present, a 
two-step process is required. One further point to note is that an equation 
analogous to (11) can be derived to relate ( z) to (^(z). 


r 
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.u»lte the above comments more precise* we will present the eqtiatiop* 
.•It if is the free space operator, i.e., Q^(z) " 0. In this 


4, (YiZ) . / I 

T 1 ‘ ir y 


(23) 


W 1 (y) = 1 

,e have 

„ / 1 sin Y z + 

7T Y 

!- r n*<n.OU - w(n)} - 3 )z - — a n 

o n - y *1 "* y 


( 29 ) 


( 30 ) 


*> 4 / ^ sin rpt + r) $(n,z) cos r,z) (l - W(n))ndn (31) 

-2 

V(r[) ~ A(n) . It is important to note that the integrands in (30) 
^i) can be shown to be even functions of r| on the real line. There- 
tne integrals may be taken from -» to » and the result halfed. Once 
done, the contour can be closed and the Residue Theorem can again 

I 

- to evaluate the integral. 

• - before, the conditions under which (27) is valid imply f(o,n) is 
/tic function of n in the upper half plane. Hence, it is sufficient 
Yfy either the magnitude or the phase. 

'*■ problem of remotely probing a stratified medium has been con- 
One formulation of the problem lias been investigated and some 
■'* results obtained. A second formulation is also presented. 

‘Jet both of these techniques are valid at arbitrary frequency 
r fi only as a parameter in the formulation), probing at frequencies 
"face waves are present requires added computation in order to 


problem. 
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However, oscillatory errors do occur in the region near z - 0 for 
many cases. It should be noted, though, that the surface fields generated 
by these erroreous profiles do not differ significantly from those generated 
by the true profile. Thus, some of the difficulty is the result of the 
poor conditioning of this type of inverse problem. 

An important feature of this method is that the numerical formulation 
requires a modest amount of computer storage so that these techniques 
might be suitable for use with measuring set-ups which have a small 
computer. This could be very helpful for diagnostic work in physics, 
biomedical applications, pollution studies and agricultural surveys. 
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Figure 2b. Inversion obtained for another linear profile. 
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, Figure 2f. Inversion obtained for a discontinuous profile. *""* 


